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Abstract 

We present an exploratory study of structure of 13 C, 15 C, 17 C and 19 C, showing that the simple 
one-quasiparticle projected BCS (PBCS) model is capable to account for several important pro¬ 
perties of these nuclei. Next we discuss the importance of the Pauli Principle in the particle-core 
models of normal-parity states in 13 C and 11 Be. This is done by considering the pairing interaction 
between nucleons moving in an over-all deformed potential. To assess the importance of pairing 
correlations in these light nuclei we use both the simple BCS and the PBCS approximations. We 
show that the Pauli Principle plays a crucial role in the parity inversion in 11 Be. It is also found 
that the effect of the particle number conservation in relatively light and/or exotic nuclei is quite 
significant. Comparison of our results with several recent papers on the same subject, as well as 
with some experimental data, is presented. 
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I. INTRODUCTION 


It has been known for a long time that the structure of the nucleus depends significantly 
on their superfluid nature. In fact, pairing constitutes the main part of the residual in¬ 


teraction beyond the Hartree-Fock 
low-energy properties of the system 


]HF) approximation and has a strong influence on most 
l|. This encompasses masses, separation energies, defor¬ 
mation, individual excitation spectra and collective excitation modes such as rotations and 
vibrations. The role of pairing correlations is particularly emphasized when going toward 
the neutron drip-line because of the proximity of the Fermi surface to the single-particle 
continuum. Indeed, the scattering of virtual pairs into the continuum gives rise to a variety 
of new phenomena in ground and excited states of nuclei Q. 

The above suggests that the pairing correlations might play a relevant role in the structure 
of light exotic nuclei. However, the fact that the n Be ground state is a l/2 + state and not a 
1/2“ state (as dictates the spherical shell model) is frequently attributed to the quadrupole 
core excitation effects only. That is, the parity inversion in 11 Be is usually described within 
the simple particle-vibration coupling model (PVM) 0, or within the simple particle- 

rotor coupling model (PRM) 7, 8|, [9j, ll(| . where the correlations among the valence particles, 
coming from the residual interaction, are totally neglected 1 . 

We know, however, that, while the positive parity states in n Be can be accounted for 
fairly well in the weak coupling model the low-lying negative parity states cannot be 
described reliably with the same model. As stated by Esbensen et al 13j, this is because 
the last have a complicated structure due to the importance of the pairing correlations. The 
same statement is also valid for the 13 C nucleus. In fact, it has been recognized since the 
paper of Lane [hlj] on 13 C in the mid-50’s that the weak-coupling of a sd nucleon to a p-shell 
core is a good starting point for the description of low-lying non-normal-parity states in a 
number of p-shell nuclei. Moreover, Lawson and Kuratli have also pointed out long ago [ul 
that the same is not true for the normal-parity states: ’’The reason is that the single-nucleon 
function which is added to the 12 C core has strong components in common with some of 


1 In the early 1960’s Talrni and Unna 11] first noticed the parity inversion in n Be and suggested that the 
interaction of a lpi /2 or 2s 1( /2 neutron with two lp 3 / 2 protons was the main cause. Many large-basis 
shell-models calculations have been carried out for n Be afterwards, but they do not all reproduce the 
parity inversion. Moreover, a recently done large basis ab initio shell model study [l3| suggests that a 
realistic 3N force will have an important influence on this regard. 
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the functions within the core.” Thus, because of the violation of the Pauli Principle (PP) 
the particle-core weak-coupling models have been mostly limited to the study the positive 
parity spectrum of 13 C j?, 16, 17]. Although important efforts have been invested to 

incorporate the PP into the core-particle models jf], 6, 18], the shell-model was so far the 
only plausible alternative to deal with the negative parity states in odd-mass carbon and 
11 Be nuclei 


19 


20 ]. 


But the pairing correlations are closely related to the PP, which is not violated within the 
BCS approximation, the fact that is not always well perceived. The problem that appears 
now, however, is the well known particle-number non-conservation. This, in turn, can be cir¬ 
cumvented through the particle-number projection procedure 


22j ], i.e., in the framework 


21 


of the projected BCS (PBCS) approximation. In fact, recently it has been found 23], 24] 
that the low-lying energy spectra in 13 C can be described quite well in the context of the 
PBCS approximation for the pairing force among the valence neutrons. Moreover, in this 
study it was shown that the projection procedure is very important to account for the weak 
decay observables around 12 C. 

The aim of the present work is twofold. First, we inquire to which extent the main feature 
of heavy carbon isotopes 15 C, 17 C, and 19 C can be interpreted within the framework of the 
PBCS model, using as a building block the results obtained in Ref. 23J, 24] for the stable 
carbon isotope 13 C. In doing so, one should keep in mind that, as the number of neutrons 
increases, in going from 13 C to 19 C, the nuclei become more and more weakly bound. There 
is no systematic study of the energy spectra of these nuclei in the literature. Secondly, we 


discuss the interplay between the single-particle and col 
n Be, by incorporating the PP into the previous works 


ective degrees of freedom in 13 C and 


correlations 


251, 126], 


4 7 

7 ^7 1 7 


10] through the pairing 


27, [28]. Therefore, we will discuss the structure of these two nuclei 


in the quasiparticle-rotor model (QPRM) and the quasiparticle-vibrator model (QPVM). 
The corresponding models with the particle number projection included will be labelled, 
respectively, as PQPRM and PQPVM. Special emphasis is put on the odd-parity states 
which so far have been treated very rudimentary within the particle-core models. 

The theoretical tools are certainly oversimplified as compared to the physical reality of 
nuclei that are discussed. In particular, neither the weak-binding characteristics nor coupling 
to the continuum states is taken into account. However, in spite of the simplicity of the 
approach it allows for obtaining simple estimates of the main structural features of these 
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nuclei. Certainly, this kind of calculation cannot replace any full-fledged description, which 
is, however, presently not available. The present analysis has to be understood as a very 
simple version of a shell-model approach, which nevertheless is able to give main components 
of the involved configurations. We also feel that the paper makes some significant advance 
in the line of research that has been quite importantly developed by previous similar studies. 


II. PAIRING DESCRIPTION OF ODD-MASS CARBON ISOTOPES 

The definitions of particle and hole states in the BCS and PBCS approximations 
are listed in Table Q] where 

e j = ifj + ^) 1/2 ; e.j = ej - Hj - l , ( 2 . 1 ) 

are the usual BCS quasiparticle energies, which depend on the single-particle energies (s.p.e.) 
ej, on the self-energy 

H = E (2j + 1)1/2 v l E Wj'\ 0 ), ( 2 . 2 ) 

and on the pairing gap 

= ~\E (2j + ijv2 u r v M:in':i'\ o). ( 2 . 3 ) 


TABLE I: Definitions of quasiparticle (Ilj + ' ) ), and quasihole (Ej ^) energies in the BCS and PBCS 
approximations. In both cases, Ej E can be either negative or positive, while Ej ^ are always 
negative. The BCS quasiparticle energies Ej, defined in (12.11) . are positive, and the corresponding 
chemical potential A is negative. The PBCS quasiparticle energies e^ are defined in (12.41) . 


Model 

d +) 


BCS 

A + Ej 

X-Ej 

PBCS 

e N 

N—2 

E 3 


4 










The PBCS energies read 




£ j = 


I N (j ) 


tdN 

W) 

IN' 


(2.4) 


The quantities R N and I N , where N is the neutron number, are defined in Ref. 21] 

The BCS and PBCS predictions for the spectroscopic factors will be also discussed 
Within the PBCS approximation these quantities are given by 

J N U) 
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S u (j) = 


(2.5) 


"3 JN ’ 

for the stripping reactions on even targets and for pick up reactions on odd targets, and by 

.I N - 2 U) 


Sv{j) = (2j + l)v i 


IN 


( 2 . 6 ) 


for the stripping reactions on odd targets and for pick up reactions on even targets. The 
plain BCS results are recovered by making all the /-factors equal to unity. 

The BCS and PBCS calculations presented here were performed in the same way as in 


the previous works 


23 


24], That is, for the residual interaction we adopted the delta force, 


V = —4Ttv pair 8(r) MeV-fm 3 , 


(2.7) 


and the configuration space includes the neutron orbital with j = nlj = 
(lsi/ 2 , 1 ^ 3 / 2 , lpi/ 2 , ld 5 / 2 , 2 si/ 2 , ld 3 / 2 ,I/ 7 / 2 , 1 / 5 / 2 )- The radial wave functions were approx¬ 
imated by that of the harmonic oscillator (HO) with the length parameter b = 1.67 fin, 
which corresponds to the estimate hto = (45/R 1 / 3 — 25A~ 2 / 3 ) MeV for the oscillator energy. 

The s.p.e. ej, as well as the value of the singlet pairing strength v pair , were fixed by 
adjusting the experimental energies, taken from 32, Table 13.4], to the calculated ones, 
through a y 2 search assuming that: a) the ground state 3/2 in n C is the lpz /2 quasihole- 


state with energy E[ p ^ , and b) the lowest 1/2 ,5/2 + ,l/2’ l ",3/2 + ,7/2 , and 3/2 states 
13 C are the quasiparticle-states j = lpi/ 2 , 1^5/2, 2si/ 2 , 1 d 3 / 2 , 1 / 7 / 2 , 2p 3 / 2 ) with energies E 


m 


(+) 


The quasiparticle energies of the distant orbitals lsi/ 2 , 2 py 2 and I/ 5 / 2 , not known 


experimentally, were assumed to be: i) E[ s ^ MeV proposed in the study done by Gillet 


and Vinli Mau 30], ii) E^ i/2 = E^ /2 + A is (p) MeV, with A is (p) = 5.73 MeV being the 
spin-orbit splitting between the p-subshclls, and iii) E^ = E^ /2 + Aho = 16.78 MeV, 
were A HO = 3.82 MeV is the energy difference between the I/ 5/2 and 2p 3 / 2 states for the 

31]. Additionally, during the minimization 


HO potential with standard parametrization 
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TABLE II: The quasihole energies Ej ^ with j = lsi/ 2 , IP 3/2 and the quasiparticle energies 
with j = lpi/ 2 , • ■ ■, I/ 5 / 2 ) used in the fitting procedure, and the resulting single-particle energies 
e,j, and the pairing strength v pair , within the BCS and PBCS for N = 6. The energies are given in 
units of MeV, and v pair is dimensionless. 


Shell 

Ef 

ej(BCS) 

ej(PBCS) 

l®l/2 

-35.00 

-23.29 

-22.27 

lp.3/2 

-18.72 

-7.74 

-7.17 

lPl/2 

-4.94 

-2.01 

-1.44 

1^5/2 

-1.09 

2.18 

2.22 

2«l/2 

-1.85 

2.78 

2.74 

143/2 

3.26 

6.88 

6.89 

1 / 7/2 

8.63 

11.04 

11.06 

2p 3 /2 

7.24 

11.61 

11.62 

2pi /2 

12.97 

17.43 

17.45 

1 / 5/2 

16.78 

19.27 

19.31 

q .pair 

u s 


23.46 

24.22 


procedure the energy difference ei P3/2 — ei Pl/2 has been kept constant and equal to A i s (p). 
The results are displayed in Table [Til 

The relationship between the s.p.e. e 3 and the corresponding quasiparticle energies 
are illustrated in Fig. [0 The major difference between the BCS and the PBCS results 
appears in and Fq P3 , which have not been used in the fitting procedure. Of course, 

opposite happens with the quasihole energies Ej ^ (see Tables III and IV in Ref. 24]). 

After fixing the parameterizations (e 3 and v pair ) in 13 C, we evaluate the low-lying energy 
spectra for the remaining odd-mass carbon isotopes, by changing only the number of neu¬ 
trons. In other words, we solve the BCS and PBCS gap equations with the parameters listed 
in Table [TTl and the number of particles N = 8,10, and 12 for 15 C, 17 C and 19 C, respectively. 
The results for the energy spectra and the spectroscopic factors are shown and compared 


with experimental data 


32 


33, 


341, 


35, 


36 


32, 


38 


39, 


40 


41 


42 


43| in Figs. [2] and [3l 
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FIG. 1: Relationship between the single-particle and the quasiparticle excitation energies for 13 C. The 
states are ordered as lsi/ 2 , lp 3 / 2 , lpi/2, 2 si/ 2 , 1^3/2, I/ 7 / 2 , 2 p 3 / 2 , 2 p x / 2 , and l/s/ 2 , and the energies are 
indicated by filled circles (BCS) and unfilled squares (PBCS). 

respectively. One immediately sees that the pairing interaction reproduces to a significant 
extent the systematic of the energy spectra of these nuclei. I 11 particular, the decrease in 
the separation energies in going from 13 C to 19 C is fairly well accounted for 2 . 

The experimentally observed spin ordering l/2 + , 5/2 + , 1/2“ of the lowest three states 
in 15 C, is well reproduced within the PBCS, indicating that the number projection plays a 
significant role, and that these levels could be predominantly one-quasiparticle (lqp) states. 
The present model also accounts fairly well for the 3/2“ state at 4.66 MeV and the 3/2 + 
state at 4.78 MeV. However, the first one of these two states could be quite likely the 
partner of the 5/2j~ state at 4.22 MeV, possessing therefore large, if not dominant, three- 

2 It would be very interesting to compare our results with a systematic large scale shell-model calculations 
of all four odd-mass carbon isotopes, employing the same Hamiltonian and the same single-particle space. 
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FIG. 2: Comparison between the calculated BCS, PBCS and measured level schemes for odd mass 


car 


m 


oon isotopes: a) BCS, b) PBCS, and c) experiments: from Ref. 


35( | for i5 C, from Refs. 


37] for 17 C, and from Ref. 


32| for 13 C, from Refs. [32|, 


33, 


36] for 19 C. For 13 C are only shown the 


experimental levels that have been used in the fitting procedure. 


quasiparticle (3qp) components | (1g? 5 / 2 ) 2 , 2; lpi/ 2 ) and |(2si/ 2 ) 2 ; 2, lpi/ 2 ), as suggested by the 


work of Bolilcu et al 


34|. 


In 15 C a doublet is known from the two-neutron transfer 13 C(t,p) 15 C reaction 44] with 
possible spin assignment of (9/2“,7/2“), at almost the same energies (6.84 MeV, 7.39 MeV) 
as the (3“, 2“) doublet in 14 C. Based on this fact, one can speculate that their 3qp structure 
is dominated by the |(1 c 7 5 / 2 ) 2 , 4; lpi/ 2 ) configuration 34]- The 1/2“ level at 5.87 MeV could 
be a 3qp (seniority one or three) state which is not contained in our configuration space. 
In fact, our restriction to lqp subspace eliminates many of the well known states in 15 C. 
However, except for the 5/2f state, the energy spectrum compares fairly well with the 
experimental results and with a recent shell-model calculation (cf. Fig. 2 in Ref. 20]). 
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The inversion of the 5/2^ and 1/2^ states in 15 C, regarding the standard single-particle 
ordering of levels, which occurs in l7 0, has been discussed a long time ago by Talmi and 
Unna [ll|. They have shown, in a brilliant manner, that the crossing of these two levels, 
in going from 17 0 to 15 C, comes from difference in interaction energy of the ld 5 / 2 and 2 si/ 2 
neutrons with two removed lpi/ 2 protons 3 . Nevertheless, as this phenomenon is typical of 
weakly-bound light neutron rich nuclei, it is seldom considered to be an “exotic” feature of 
15 C |35l. 145]. That is, it is attributed to the halo formation, which makes the lowest angular 
momentum to gain energy by extending its wave function. This apparent “anomaly” even 
takes place in the non-exotic 13 C nucleus, and in our model it can be interpreted as a 
consequence of the pairing interaction. In fact, as seen from Table HT1 and Fig. El we obtain 
that for all odd carbon isotopes is > E^ /2 , although e±d 5/2 < e 2si/2 . This is a direct 

consequence of fact that the 2 si/ 2 self-energy is significantly larger than that of the lc? 5 / 2 
state, because of the strong interaction between the lsi/ 2 and 2 si/ 2 states (see Eq. (12.21) 1. 
It might be worthwhile to point out that the effective s.p.e. e-j used here should not be 
confounded with the 17 0 energy spectrum. 

The separation energy of the last neutron in 17 C is S n = 729 ± 18 keV 46], and the shell 


model calculation 


47] predicts a J n = 3/2+ ground state. This prediction has been con¬ 


firmed later on by the single-neutron knockout reaction measurements done by Maddalena 


et al 


39], which strongly indicates such an assignment instead of the naively expected option 


J 71 " = 5/2 + arising from the seniority-one state |(ld 5 / 2 ) 3 J 7r = 5/2+). A simple explanation 
for this experimental result could be found in the so called J = j — 1 anomaly discussed by 
Bohr and Mottelson 48]. In fact, since the work of Kurath 49] we know that for (j) 3 con¬ 
figurations with j > 5/2, the J = j — 1 state can occur below the J = j state for sufficiently 
long range forces 4 . The model discussed here does not contain seniority-three states and 
therefore it is unable to account for the 3/2 + ground state spin. However, as seen from the 
Fig. [2], it predicts that the first two excited states are l/2 + and 5/2+ (in this order), which 


3 The fact that 17 C has a low-lying 1/2+ state, and 19 C has a 1/2+ ground state could also be a consequence 
of the same lp^/ 2 ld 5 / 2 and 1 p^/ 2 2si/ 2 (T = 1) particle-hole matrix elements, lowering configurations 
involving the 1/2+ states in 19 0 and 21 0 at 1.47 MeV and 1.22 MeV, respectively, very much in the same 
way as for 15 C. 

4 The 3/2+ ground state in 17 C can similarly be attributed to the differential effect of above mentioned 
particle-hole interaction on the 2s!/2 content of the dominantly (lds/ 2 ) 3 5/2+ and 3/2+ states of 19 0. 
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FIG. 3: (Color online) One-particle reaction spectroscopic factors as a function of the mass number 
A. Upper panel: S u (j ), for stripping on even parent and pick-up/ knock-out on odd parent. Lower 
panel: S v (j), for pick-up/knock-out on even parent and stripping on odd parent. The dashed and 
solid lines correspond, respectively, to the BCS and PBCS predictions. The experimental values for 
the lowest J n = 1/2“, 3/2“, l/2 + and 5/2 + states are indicated, respectively, with circles, triangles, 


squares, and diamonds. The data were by recompiled from Refs.: 
and 42] [el . 


3i [al , halide 


,QEk 


421 


36]. 


is consistent with the recent measurement done by Elekes et al 

Special attention was given to the neutron-rich carbon isotope 19 C to establish whether it 
has a pure S 1/2 one-neutron-halo as suggested in Refs. |1Q|, 150]. On the basis of measurements 


of different observables, associated with the neutron-removal reaction 


39 


51 


52], there 


seems to be the consensus that the spin and parity of its ground state is J n = l/2 + . 
Contrarily, there is a strong discrepancy in the literature regarding the separation energy of 
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19 


C. The tabulated values go from S n = 160 ± 95 keV in 1993-1997 53] to S n = 580 ± 90 


keV in 2003 


54] . The experiments using time-of-flight techniques suggest small separation 


energy, that is, weighted average yiclc 


s S n — 242 ± 95 keV 50]. The Coulomb dissociation 


of 19 C was studied by Nakamura in 51], and the analysis of angular distributions of breakup 


products suggests the value S n = 0.53 ± 0.13 MeV. Using this value in the simple cluster 
model calculation of the dipole strength gives good a gree ment with the data. On the other 


hand, more recent experiment of Maddalena et al. 


39] on nuclear breakup of 19 C yields 


S n = 0.65 ± 0.15 MeV and 0.8 ± 0.3 MeV. In the present calculation we correctly reproduce 
the spin and parity of the ground state. For its energy we obtain —0.33 MeV in the BCS 
and 0.12 MeV in the PBCS. Both results are in fair agreement with the values reported in 


Refs. 


50 


53]. Two excited states at energies of 197(6) keV and 269(8) keV were reported 


in a recent study of the y-ray spectra 


36]. The suggested spins and parities are 3/2 + 


and 5/2 + , and, as seen from Fig. [2j we reproduce the second one only. The level 3/2 + is 
very likely dominated by the seniority-three configuration |(ld 5 / 2 ) 3 </ 7r = 3/2 + ), same as the 
ground state in 17 C. 

The BCS and PBCS one-particle reaction spectroscopic factors for the lowest J 77 = 
1/2 - , 3/2~, l/2 + and 5/2 + states in odd-mass carbon isotopes are shown in Fig. [3] as a 
function of the mass number A. They are quite similar to each other, and, except for 13 C, 
they also agree fairly well with the experimental data which are displayed in the same figure. 
This agreement clearly implies that in 15 C- 19 C nuclei the states J 77 = 1/2 - , l/2 + and 5/2 + 
are basically seniority-one states. In contrast, the contribution of seniority-three states, 
and/or of the collective degrees of freedom seems to be quite relevant in 13 C. The latter will 
be discussed in the next section. 


III. i3 C AND n BE WITHIN THE QUASIPARTICLE-CORE MODEL 


S,E 


The quasiparticle-rotor model hamiltonian is evaluated in the same way as in the PRM 


a 


lot, except for: 


1. The band-head energies are modified as: 


e i E 0 


(+) 


(3.1) 
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with E ( j +> defined in Table [Q and 

2. The non-diagonal particle-core matrix elements are renormalised by the overlap factors: 

F ?y - u i u y ~ v i v E ( 3 - 2 ) 


in the QPRM, and 


in the PQPRM. 


rp p Q p _ 


V (.//) - W /- v 

[I N (j)I N U ')] 1/2 


(3.3) 


TABLE III: Neutron Wood-Saxon s.p.e. ej for 13 C and 11 Be and the corresponding quasiparticle 
BCS and PBCS energies e[ F . together the energies E(2 + ) of the collective 2 + states. All energies 
are given in units of MeV, and the pairing strength is v^ air = 30. 


Shell 

13 C 

n Be 


e 3 

BCS 

Ej + » 

PBCS 

e i 

BCS 

E y> 

PBCS 

1^3/2 

-13.54 

-9.5956 

-7.109 

-7.33 

-3.5928 

-1.8592 

Rh/2 

-7.82 

-10.4057 

-12.209 

-2.59 

-3.9583 

-5.4577 

1^5/2 

-0.82 

-4.2699 

-4.123 

3.50 

0.2640 

0.2333 

2-S1/2 

-0.61 

-3.4033 

-3.303 

0.30 

-2.0479 

-2.1316 

E{ 2+) 

4.438 

3.368 


Thus, the PRM differs in several important aspects from the QPRM and PQPRM. First, 
as shown in Fig. [Q the BCS and PBCS energies Ej^ can be quite different from the s.p.e. 
ej , not only in magnitude but also in sign. Second, the factors Fjj>, which are quite similar 
in BCS and PBCS, correctly take in account the PP and they can by considerable less than 
unity for states near the Fermi level, diminishing in this way their coupling to the core 
excited states quite a lot. In addition, the particle-like states do not couple to the hole-like 
states, and if the particle-core coupling is attractive (repulsive) for a particle-like state j, it 
is repulsive (attractive) for a hole-like state with same quantum numbers j. 

The vibrational analog of the QPRM, i.e., the quasiparticle-vibrator coupling model 
(QPVM) has been introduced in Refs. [27, 28, 4s|. Note also that when only the collective 
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FIG. 4: (Color online) Calculated levels for 13 C as functions of the /3 value for 12 C, within the 
PRM, QPRM and PQPRM models. 


state 2 + is considered, and the diagonal 2 + — 2 + interactions are neglected, the QPRM and 
QPVM formally yield the same result, except that the negative value of the deformation 
parameter (3 does not have any physical meaning in the vibrational case. The same statement 
is valid for the projected version of the QPVM, i.e., for the PQPVM. 

Incidentally the s.p.e., used in the previous section and shown in Table UT1 were extracted 
from the experimental data and thus already include the collective degrees of freedom. As 
such, they can not be used in the core-particle models. Moreover, due to the fact that the 
results strongly depend on the size of the configuration space, and in order to make the 


comparison with previous particle-core calculations 


,14 


7, 


10] as close as possible, 


we will use here only four bare s.p.e., obtained from a Wood-Saxon potential with standard 
parametrization 55]. The results for 13 C and n Be are shown in Table Hill Being the single- 
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PQPRM EXP PVM PRM 


FIG. 5: (Color online) The low-energy spectra in 13 C within the PQPRM for (3 = —0.6, compared 
with the experimental levels (EXP) and the PVM [3] and PRM [{| calculations. 


particle space here smaller than the one used in the previous section, we employ now a 
somewhat larger pairing strength, in order to achieve the convergence of the BCS equations. 
We wish to stress that all the discussion done within the particle-core models is basically 
qualitative, and therefore the choice of the model parameters (Wood-Saxon potential, E( 2 + ), 
uf a * r , and (3) doesn’t play a crucial role here. 

The calculated low-lying levels for 13 C as functions of the (3 value for 12 C, within the 
PRM, QPRM and PQPRM models are shown in Fig. [4j Let us first note that the lowest 
states P = 3/2R, 1/2R, 3/2^ and 5/2^, based on the unperturbed configurations |lp 3 / 2 , 0+) 
and |lp 3 / 2 ,2 + ), are frequently just ignored within the particle-core coupling models [3j, 4, 
8, sl in order to simulate the PP 5 . We are considering them, however, in order to make 
the comparison with the QPRM and PQPRM models, where these states are physically 
meaningful. Moreover, the state P = 1/2^, build up on |lpi/2,0 + ), is considered to be the 


5 In Ref. [7J have not been considered the negative parity states. 
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ground state in the PRM. 

The measured quadrupole moment of the core 12 C is Q 0 = — (22±10) e fm 2 56|, suggesting 
an oblate shape. This gives a quite large quadrupole deformation (3 = —0.6). As can be seen 
from Fig. SJ only the PQPRM reproduces satisfactorily the experimental energy ordering of 
the lowest four levels in 13 C with this value of 3. In Fig. Owe confront the PQPRM energies 
for this 3 °f the lowest l/2~, 3/2“, l/2 + and 5/2 + levels in 13 C with the experimental data, 


and the PVM [3| and PRM 9] calculations. 

It should be noted that the PRM and QPRM are unable to account for the 13 C energy 
spectra for any value of the deformation parameter, neither positive nor negative. In the 
first case, and even when the levels based on the configurations |lp 3 / 2 ,0 + ) and |lp 3 / 2 ,2 + ) 
are omitted, a low lying 5/2~ state, arising from the unperturbed |lpi/ 2 , 2 + ) level, shows up, 
which is not observed experimentally. On the other hand, in the second case, the particle- 
core coupling never removes the degeneracy between the states 1/2R and 3/2R. 

The resulting PRM, QPRM and PQPRM wave functions for the lowest 1/2A 3/2" 
l/2 + , and 5/2 + states, with 3 = —0.6, are confronted with previous calculations 


BE, 


m 


Table HV1 As mentioned before, within the PRM are ignored the states J 77 = 3/2R, 1/2R, and 
3/27, based on the unperturbed configurations |lp 3 / 2 ,0 + ) and |lp 3 / 2 ,2 + ), and are exhibited 
the wave functions of the l/2 2 and 3/2^" levels. On the other hand, in the PRM calculations 
done in Refs 7, 9] instead of the wave functions are presented the percentages of the different 
configurations, and therefore we only show absolute values of the corresponding amplitudes. 

It is worth noting that, while the pairing makes the wave functions of the negative parity 
states less collective, all three models used here yield similar results for the positive parity 
states. We find that about 55% of the 1 /2^" state consists of the ld 5 / 2 single-particle state 
coupled to the 2 + excited core state. In contrast, in the previous particle-core coupling 
calculations |3|, 7, 9] this state is basically (> 90%) of single-particle nature. As seen from 
the Table [IVl our wave functions for the 5/2|" are also significantly more collective that 
those obtained in the just mentioned works. 

The differences in the wave functions, displayed in the Table IIV1 are reflected in the 
reaction spectroscopic factors for pickup on the 13 C target and for stripping on the 12 C 
target, which are shown in the same table. Also are listed the results of the previous PVM 


and PRM 


BBS, 


and shell model [19( studies, as well as the experimental results 


42 


3 , 


57J; 


which are accounted for quite well within the present PQPRM calculations. 
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TABLE IV: Comparison between the wave functions and the spectroscopic factors obtained in the present 
work for the lowest states in 13 C with those derived by Vinh Mau 3] within the PVM, and by Nunes et al 

q n ' 

9j and Esbensen et al [7J within the PRM, and the shell model study done by Cohen and Kurath 
all our calculations (PRM, QPRM, and PQPRM) we adopt (3 = —0.60. 


19J. In 



Model 

Wave Functions 

Spectroscopic Factors 



|ipi/2;0 + >|ipi/ 2 ;2+> |ip 3 / 2 ;0 + ) |ip 3/2 ;2 + > 

Theor 


Exp. 








Ref. [57] 

Ref. [43] Ref. [42] 

1/2- 

PRM 

0.849 



0.528 

0.72 




QPRM 

0.997 



0.076 

0.69 




PQPRM 

0.961 



0.277 

0.53 

0.58 (15) 

0.77 0.61 (9) 


PVM [3] 

0.791 



0.602 

0.63 




PRM [9] 

0.565 



0.818 

0.32 




SM \W] 





0.61 



3/2" 

PRM 


0.929 

-0.162 

0.332 

0.03 




QPRM 


-0.052 

0.935 

0.350 

0.17 




PQPRM 


0.759 

-0.650 

0.025 

0.10 


0.14 


PRM [9] 


|0.787| 

0.375 

0.489 

0.14 




SM \W] 





0.19 





|2si/ 2 ; 0+) |2s!/ 2 ; 2+) |ld 5 / 2 ;0+) ld 5 / 2 ;2+) 


1/2+ 

PRM 

0.684 



0.730 

0.47 




QPRM 

0.672 



0.741 

0.44 




PQPRM 

0.671 



0.741 

0.44 

0.36(2) 

0.65 


PVM [3] 

0.957 



0.291 

0.92 




PRM [7] 

0.975 



0.200 

0.95 




PRM [9] 

0.952 



0.272 

0.90 




SM [58] 

0.945 



0.251 

0.89 



5/2+ 

PRM 


0.485 

0.737 

-0.471 

0.54 




QPRM 


0.463 

0.752 

-0.469 

0.55 




PQPRM 


0.466 

0.749 

-0.471 

0.55 


0.58 


PVM [3] 


0.179 

0.867 

0.300 

0.75 




PRM [7] 


0.144 

0.894 

0.424 

0.80 




PRM [9] 


0.565 

0.831 

0.538 

0.69 




SM [58] 


? 

0.897 

-0.357 

0.80 




We obtain much smaller amplitudes for the non-collective components |2si/ 2 ; 0 + ) and 
|ld 5 / 2 ;0 + ) of the lowest |l/2[/) and 15/2i”) states than do the full 1 hco shell-model calcula¬ 
tions, among which the work of Jager et al 59] was very likely the first one. There have 
been many such calculations since then with refined effective interactions, but the wave 
functions for the lowest positive-parity states remain essentially unchanged. We show the 
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ones obtained in the work done by Millener et al 58j, which, as seen from Table IIVI yield 
spectroscopic factors that are ~ 30 — 60% larger than the measured ones. 

In Fig. Owe show the PRM, QPRM and PQPRM results of our study of the low-lying 
states in 11 Be as a function of deformation. All what has been pointed out in commenting 
the 13 C nucleus in the PRM, regarding the PP, the unperturbed configurations |lp3/2,0 + ) 
and |lp 3 /2,2 + ), and the negative parity states, is also pertinent here. On the other hand, 


0 


when comparing our PRM results with the work of Esbensen et al [7H one sees that their 
deformation dependence of the levels l/2 + and 5/2 + gets flatter than ours for strong defor¬ 
mations. The reason for that is simple and comes from the fact that they use a volume factor 


in the calculation, given by [7, (3)], to preserve the volume of the nucleus. This actually 
means that their interaction radius is smaller for higher deformations and this results in the 


d be noted that our PRM curves for the 
, Fig. 2], The main difference is that in 


flattening of the curves. Having said that, it shou 
states 1/2/ and 1/2^ are very much like those in 
our case the crossings between the positive and negative 1/2 states occur at a significantly 
smaller value of 3. 

Both PRM calculations performed so far [7, 9] were done with a positive value of f3, i.e., 
a prolate deformation has been assumed for the 10 Be nucleus. However, there is no firm 
experimental evidence that it is so. Furthermore, we neither know whether 10 Be behaves as 
a rotator or as a vibrator. From its energy spectra, with the first and second 2 + at energies 

= 3.37 MeV and E~+ = 5.96 MeV (i.e., E~+ = 2 E 9 +), one can conclude that it is more 
i 2 n 2 1 

likely a vibrator. In fact, Esbensen et al [7] have noted that for the positive parity spectra in 
11 Be a better agreement with data is obtained after reducing (in 34%) the 2 + — 2 + coupling 
strength. On the other hand, Vinh Mau {s| and Colo et al {4 ] use straightforwardly the 
PVM. 

The most recent study of quadrupole deformation of 10 Be has been done through proton 


inelastic scattering 60], and the values for /3, extracted from the measured deformation 
lengths 5 = PR = Pr^A 1 ^, are: P = 0.593(56) and P = 0.692(65). A careful scrutiny of 
the energy levels in Fig. [6l leads one to conclude that none of our particle-rotor coupling 
calculations is able to reproduce the experimentally observed spin sequence l/2 + — 1/2“ — 
5/2 + in n Be, for such a value of p neither positive nor negative. In view of this it does not 
make much sense to comment the wave functions within the PRM, QPRM and PQPRM 
approximations. 
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FIG. 6: (Color online) Calculated levels for 11 Be as functions of the (5 value for 10 Be, within the 
PRM, QPRM and PQPRM models. 

The just mentioned results have induced us to replace the rotor by a harmonic vibrator in 
the description of the low-lying n Be levels. The results are shown in Fig. [7] The PVM (upper 
panel) exhibits the same features as the PRM regarding the PP. Therefore the lowest four 
negative parity states, same as in the PRM, have to be discarded within the PVM. However, 
even doing so it is not possible to obtain the desired theoretical results. On the other hand, 
from the QPVM energy spectra (middle panel) one sees that the particle-vibrator coupling, 
same as the particle-rotor coupling, is unable to remove the degeneracy between the states 
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p 

FIG. 7: (Color online) Calculated levels for 11 Be as functions of the (5 value for 10 Be, within the 
PVM, QPVM and PQPVM models. 


1/2 and 3/2 jb Once more, as seen from the PQPVM results (lower panel), this is achieved 
through the number projection procedure only. 

Both the 1/2/ and 1/2 j" levels go down within the PQPVM when /3 is increased, and their 
crossing happens close to the experimental value for the deformation parameter (j3 = 0.6). 
Thus, the spin inversion mechanism is quite different here than in the PRM and PVM cases, 
where the 1/2state is pushed up while the 1/2/ state is pushed down jl], s|. The fact that 


in the PQPVM the 1/2 1 state is almost independent on the value of f3 while the 1/2/ state 
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3 - 


2 




1 - 


3/2 
5/2 + 


1 / 2 - 


1/2-i 


PQPVM 


( 1 / 2 , 3 / 2 , 5 / 2 + ) 111 


( 1 / 2 , 3 / 2 , 5/2 ) 



. 3 / 2 - 


5/2 + 
— l/2~ 


FIG. 8: (Color online) The low-energy spectra in n Be within the PQPVM for (3 = 0.6, compared 
with the experimental levels (EXP) and the PVM [3] and PRM [fj| calculations. 


varies very strongly is due to the PP factor given in the equation (13.31) . which is small in 
the first case (u Pl/2 = v Pl/2 ) and large in second case {u Pl/2 = l,u Pl/2 = 0). Therefore, it can 
be argued that within the particle-core model the spin inversion in 11 Be is due to the PP. 

As shown in Fig. [8] all three lowest states in 11 Be are satisfactorily reproduced within the 
PQPVM for (3 = 0.6. For the sake of comparison in the same figure are also presented the 
results of the PVM and PRM 9] calculations. 

The resulting PVM, QPVM and PQPVM wave functions for the lowest 1/2 - , 3/2 _ , 1/2 + , 
and 5/2 + states in 11 B, with f3 = 0.6, are confronted with previous calculations 3, 7], |9, j61] 
in Table IVl All said above in relation of omission of levels based on configurations |lp 3 / 2 , 0 + ) 
and |lp 3 / 2 , 2 + ), and the percentages of the wave function amplitudes in the case of 13 C is 
valid here. 

The PVM wave function for the 1/2R state, obtained by Vinh Mau jsj], is quite similar 
to the one we get within the same model. Notice that the pairing tends to make this state 
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TABLE V: Comparison between the wave functions and the spectroscopic factors obtained in the present 

TinfVi fli aoq /^nvnriarl Vvit - AfvnVi A/Tfin a withir 


work for the lowest states in 11 Be with those derived by Vinh Mau [3j within the PVM, and by Nunes et 


a 


al 9fl and Esbensen et al [7} within the PRM, the variation shell model (VSM) of Otsuka et al 
the shell model study done by Cohen and Kurath 
adopt P = 0.60. 


61j . and 


19|. In all our calculations (PVM, QPVM, PQPVM) we 


J n 

Model 

Wave Functions 

Spectroscopic Factors 



|ipi/2;0 + )|ipi/ 2 ;2+)|ip 3 / 

2 ;0+) |lp 3 / 2 ;2 + ) 

Theor 


Exp. 









Ref. [62] 

Ref. [63] Ref. [64] 

1/2- 

PVM 

0.750 



-0.631 

0.56 





QPVM 

0.995 



-0.103 

0.63 





PQPVM 

0.916 



-0.401 

0.49 

0.63 (15) 

0.96 

0.45 (12) 


PVM [3] 

0.746 



0.667 

0.56 





PRM [9] 

0.932 



0.362 

0.87 





SM [19] 





0.60 




3/2" 

PVM 


0.870 

0.428 

-0.247 

0.18 





QPVM 


0.070 

0.927 

-0.368 

0.25 





PQPVM 


0.592 

0.735 

-0.331 

0.15 





SM [19] 





0.17 






2s 1/2 ; 0+) |2s 1/2 ; 2+) |1g? 5/ 

2 j 0+) |irf 5 / 2 ; 2+) 


1/2+ 

PVM 

0.826 



-0.564 

0.68 





QPVM 

0.825 



-0.565 

0.62 





PQPVM 

0.821 



-0.571 

0.64 

0.73 (6) 

0.77 

0.53 (13) 


PVM [3] 

0.964 



0.267 

0.93 





PRM [7] 

0.933 



0.317 

0.87 





PRM [9] 

0.883 



0.447 

0.78 





VSM [6lj 

0.740 



0.630 

0.55 




5/2+ 

PVM 


-0.523 

0.740 

0.422 

0.54 





QPVM 


-0.486 

0.760 

0.375 

0.56 





PQPVM 


-0.485 

0.758 

0.437 

0.56 


0.50 



PVM [3] 


0.269 

0.896 

0.353 

0.80 





to be less collective. The effect is still more pronounced for the 3/2 1 state. 

All three models yields similar results for the positive parity wave functions. When 


confronted with that derived within the PVM by Vinh Mau [ 2 ]], one sees that, at variance 
with what happens with negative parity states, our wave functions are more collective. 
We also note that our wave function for the l/2/~ state is similar to that obtained in the 
variational shell model (VSM) calculation of Otsuka et al Ref. 61]. 

The reaction spectroscopic factors for pickup on the n Be target and for stripping on the 
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reproduced quite well by the PQPVM. 


Before ending this section, and in order to inquire on how realistic the above particle-core 
wave functions are, we briefly discuss the ground state magnetic dipoles moments in 13 C and 
11 Be. After some Racah algebra and by denoting a and b , respectively, the single-particle 
and collective components in the wave functions, we get: 

MV 2 r) = a 2 g{lp 1/2 ) - yhOs/2) + yh( 2+ ) 

= ~^r + b 2 gR, (3.4) 

o 


and 


7b 2 b 2 

Ml/2^) = a 2 g(2s 1/2 ) + —g(ld 5/2 ) - —g(2 + ) 


y± 

2 


a- + 


15 

7 b 2 \ 2 b 2 

15 


-Qr, 


(3.5) 


where g s = —3.82 and g r = Z/A are, respectively, the effective spin and collective gyromag- 
netic ratios 481,1551]. 


TABLE VI: Magnetic dipole moments of the ground states in 13 C (1 /2j“) and n Be (I /2f ) evaluated 
with the wave functions listed in Tables HVl and |V]are confronted with the experimental data. PCM 
stands here for the particle-core model, with the core being a rotor in the case of 13 C, and a vibrator 
in the case of n Be, etc. 


Nucleus 

Theory 

Experiment 


PCM QPCMPQPCM Ref. [3] Ref.[9j Ref. [7] Ref. [61] 

Ref. [65] 

13 C 

0.77 0.63 0.67 0.80 0.94 

0.7024118(14) 

n Be 

-1.66 -1.67 -1.66 -1.85 -1.78 -1.71 -1.49 

-1.6816(8) 


The numerical results are shown in Table IVH from where one can see that they the 
magnetic dipole moments strongly depend on the interplay between the single-particle and 
collective degrees of freedom. 
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IV. SUMMARY AND CONCLUSIONS 


Our work was divided into two main stages. First, we apply pure BCS and PBCS 
models to describe heavy carbon isotopes. We adopt the single-particle energies and pairing 
strengths which reproduce simultaneously the experimental binding energies of 11 C and 13 C, 
and the low-energy spectra of 13 C, except for the 3/2f state, which is nicely reproduced 
within the PBCS but not within the BCS (see Fig. [2]). With these parameters we next 
evaluate the low-lying states in 15 C, 17 C and 19 C. We found that both models are capable 
to explain fairly well the decrease of the binding energies in going from 13 C to 19 C. The 
PBCS model reproduces as well the experimentally observed spin ordering l/2 + , 5/2 + , 3/2 + 
of the lowest three states and the energy of the first 3/2“ state in 15 C. In addition, the same 
model correctly predicts the lowest l/2 + and 5/2 + levels in 17 C and 19 C, but it is not able 
to account for the 3/2^ level in these two nuclei, which is very likely build up mainly on the 
anomalous seniority-three state |(1 d 5 / 2 ) 3 J 7T = 3/2 + ). 

We have calculated the one-particle transfer spectroscopic factors for the lowest 
l/2~, 3/2“, l/2 + and 5/2 + states for all four odd-mass C isotopes obtaining quite similar 
results within the two pairing models. In 13 C the calculations agree with the experimental 
data for the ground state, but they fail to reproduce the experimental spectroscopic factors 
for excited states in this nucleus. For the heavier carbon isotopes the agreement with the 
experimental data is better. This implies that the low-lying states in 15 ’17 ,i9q are basically 
seniority-one states, while the contribution of seniority-three states and the collective effects 
could be important in 13 C. Our results clearly indicate that the pairing interaction plays a 
major role in the nuclear structure of heavy carbon isotopes, partly accounting for their ba¬ 
sic feature such as: a) small binding energies, b) spin-parity ordering of the low-lying states, 
and c) systematic decrease in the binding (of the order of 5 MeV) when one goes from 13 C 
to 19 C. All this is achieved without employing any free parameter. Therefore, the pairing 
has to be taken into account in any theoretical calculation that aspire to be quantitatively 
realistic. 

In the following stage we included the collective degrees of freedom in the framework of 
the weak-coupling model and applied it to describe 13 C and 11 Be. Our main objective here 
was to analyze how the short range pairing correlations modify the core-particle coupling 
mechanism, and consequently the energy spectra, spectroscopic factors, and the magnetic 
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dipole moments in these nuclei. As far as we know such a study has not been done so far, at 
least not in a systematic way 6 . It is worthwhile to stress once more that the PP, which is 
usually omitted in the simple particle-rotor and particle-vibrator coupling models, is brought 
up back into the play by the inclusion of the pairing. Here the single-particle energies were 
taken from the standard parametrization for the Wood-Saxon potential, while the coupling 
matrix elements were calculated using the wave functions of the harmonic oscillator. The 
last procedure could be a rather crude approach for a weakly bound nucleus such as 11 Be. 
However, we feel that it is, nevertheless, good enough to reveal the importance of the PP 

We found that only PQPRM reproduces satisfactorily the experimental ordering of the 
lowest four levels in 13 C with the correct value of deformation (3 = —0.6. Neither the PRM 
nor the QPRM could accomplish this. The pairing strongly reduces the coupling between the 
0 + and the 2 + core states for the negative parity states 1/2“ and 3/2“. As a consequence, 
within the QPRM and PQPRM the ground state of 13 C turns out to be basically the single¬ 
particle lpi /2 state. Our positive parity states l/2 + and 5/2 + , on the other hand, are more 
collective when compared to results of the previous works. Both effects seems to go in right 
direction and make the PQRPM to account for the experimental spectroscopic factors of 
low lying positive and negative parity states, and for the dipole magnetic moment of the 
ground state. 

Similar calculations were done for one-neutron halo nucleus 11 Be. It is found that the 
more likely structure of the core, 10 Be, required to reproduce the lowest three states in 11 Be, 
is that of a vibrator. Moreover, we feel that the appropriate model for the experimental 
value of the vibrational length, (3 = 0.6, is again the PQPVM. As before, we found that 
inclusion of pairing makes the negative-parity states less collective and positive parity states 
more collective compared to earlier works. Also here the spectroscopic factors and the dipole 
magnetic moment are well reproduced by the PQPVM. 

In summary, the inclusion of the pairing interaction and of the concomitant Pauli Prin¬ 
ciple is imperative not only in the case of heavy odd-mass carbon isotopes but also in the 
core-coupling models of 13 C and n Be. The important role played by the particle number 
conservation in these relatively light and/or exotic nuclei has been confirmed as well. 

6 See however the Refs. B lii ITsl . 
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